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Abstract 

The semiclassical Coulomb excitation interaction is at times expressed in the Lorentz gauge in 
terms of the electromagnetic fields and a contribution from the scalar electric potential. We point 
out that the potential term can make spurious contributions to excitation cross sections, especially 
when the decay of excited states is taken into account. We show that, through an appropriate 
gauge transformation, the excitation interaction can be expressed in terms of the electromagnetic 
fields alone. 
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Introduction 



Coulomb excitation has proven itself as an important tool for studying the structure of 
both stable and exotic nuclei . The phenomenon has been well-studied in a nonrel- 

ativistic context Q| and studied perturbatively in a relativistic onep] by Alder and Winther. 
More recent experimental studies have shown the need to consider multiple Coulomb exci- 
tatioa at relativistic eae rg ie S fl Q Q aad, althoa gh tae fc, m ali Sm of Kef. caa be extended 
to permit such calculations, it is not easy to implement 0, Isj]. An alternative semiclassical 
form of the Coulomb excitation interaction was presented in Ref. Is and refined in a more 
recent analysis |lo|. There, the Coulomb excitation interaction is expressed in the Lorentz 
gauge in terms of the electromagnetic fields and the scalar electric potential. This mixed 
representation of the interaction yields the results of Winther and Alder when used pertur- 
batively and is adequate for describing single and multiple excitation of states of zero width 
at incident energies lower than about 2 GeV/nucleon, where the results are almost identical 
to those obtained by directly using the electromagnetic fields, as in a classical treatment of 
the problem. At higher energies, coupled channel calculations in the mixed and field rep- 
resentations yield results that are increasingly discrepant, as the excitation cross section of 
the mixed representation quickly grows to absurdly large values. Similar discrepancies have 
also been observed in a recent comparison of relativistic Coulomb excitation in the Lorentz 
and Coulomb gauges 

Electromagnetic processes, such as Coulomb excitation, should be gauge-invariant. Trun- 
cation of the coupled-channels model space can spoil this invariance however permitting 
spurious, but gauge-removable terms, such as the potential term in the Lorentz interaction, 
to contribute to the cross sections. Baltz, Rhoades-Brown and Weneser have seen similar 
,a rg e effects ia the, extensive stnd y o, e +e - prodnctionQ U Q aad nave fonnd as well 
that, when performed with care, an appropriate gauge transformation can greatly simplify 
calculations. 

In the case of the mixed interaction of the Lorentz gauge, the effects of the spurious 
potential term are exacerbated even further when the excited states are allowed to decay. 
The mixed representation then yields large, unphysical cross sections at all incident energies, 
due to the loss of flux after excitation by the long-range potential term. The production 
cross section of the decay products of a dipole transition, in particular, grows as the square of 
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the Lorentz factor 7. This makes a treatment of multiple Coulomb excitation incorporating 
fluctuation contributions, such as those of the Brink- Axel type [3, 0, Q, Q| , unviable in 
the mixed representation, since finite widths are a fundamental component of such models. 

In the following, we first provide an estimate of the production cross section of the de- 
cay products of a dipole transition induced by the long-range potential term. We then 
demonstrate, through an appropriate gauge transformation, that the Coulomb excitation 
component of the interaction can indeed be expressed in terms of the electromagnetic fields 
alone. Aside from making a satisfying parallel with the classical case, the pure field rep- 
resentation of the excitation interaction, known as the multipolar or Poincare gauge) 20], 
provides physically reasonable cross sections even when the excited states decay[18]. A sim- 



ilar expression has 
in metallic clusters 



3een used in the treatment of Coulomb excitation of plasmon resonances 
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The Semiclassical Coulomb interaction 

In the usual semiclassical approximation to heavy-ion scattering, the relative motion is 
described by a classical trajectory. The projectile-target interaction is then a time-dependent 
function determined by this trajectory. The semiclassical form of the electromagnetic inter- 
action Hamiltonian is given by 

V(t) = J d 3 x (p(x,t)<p(x,t) - ^J(x,t) ■ A(x,t)^j , (1) 

where <f(x,t) and A(x, t) are the scalar and vector electromagnetic potentials due to the 
projectile, for which 

E(x, t) = -V<p(x, t) - 1 ^4(M) and t) = V x A(x, t), (2) 

c o t> 

and p(x, t) and J(x, t) are the charge and current density operators of the target nucleus. In 
high-energy collisions, one usually uses the Lienard-Wiechert potential, which is the retarded 
electromagnetic potential, in the Lorentz gauge, of a charged particle moving on a straight 
line trajectory. 

In Ref. 0, the Lienard-Wiechert potential was used to obtain the Coulomb interaction 
for multipole excitation of a target nucleus in a mixed representation that depends on both 
the electric field and the scalar electric potential. The interaction inducing transverse dipole 



excitations is written in terms of the electric field alone. Due to the contribution of the 
vector potential, the interaction inducing longitudinal transitions also includes a potential- 
like term. It can be written as 



Vn(b,t) = V 



■iut 



(3) 



(b 2 + (jVt) 2 f /2 C dty (&2 + ( 7 ^)2)V2 

with b being the impact parameter, v the projectile velocity, huo the excitation energy, 
(3 = v/c and 7 the associated Lorentz factor. The quantity Vq represents the product of the 
projectile charge, the matrix element for dipole excitation and numerical factors. (Note that 
the factor £z{t) defined in Eq. (26) of Ref. must be divided by r 2 to provide the correct 



expression, which can be found in Ref. 



Let us now consider the time evolution of three states: the ground state, the longitudinal 
dipole mode excited by Coulomb excitation and the residual state fed by the decay of the 
latter. The time evolution equations of the ground-state and dipole-mode amplitudes, a (b, t) 
and ax\\(b,t), respectively, can be written in the interaction picture as[3| 

ih^a (b,t) = Vu{b,t)e- itjJt a n (b,t) 
dt 

d r 

ih-a n (b,t) = ^(M^VM) - i-a n (b,t), (4) 

where T is the width of the dipole resonance. The residual state is fed incoherently by the 
decay of the dipole mode. The time evolution equation for its occupation probability can 
be written as 

%(M)^MM)f. (5) 

Coupled coherent-incoherent evolution equations such as these can be consistently formu- 
lated in terms of the density matrix, as shown in Refs. Q and Q. The formulation above 
is sufficient for our purposes here. 

Since the interaction tends to zero as b — > 00, first order perturbation theory is valid at 
large values of the impact parameter, where depletion of the ground-state and occupation 
of the longitudinal excitation can be neglected on the right side of the equation. We can 
then approximate the amplitude for excitation of the longitudinal mode by retaining only 
the second term in Vx\\(b,t), which decreases as as 

ai||(M) ~ -t / dsV 1 \\(b,s)e VJS (6) 

= IVq- -pr + 0\ — 

he (fo2 + (jvt) 2 ) 1/2 \b 3 



It is clear that the contribution of the 1/6 term vanishes as t — > oo and that it thus makes 
no net contribution to the excitation, when no flux is absorbed from the excited mode and 
ground-state depletion is negligible. These are the assumptions used in the usual perturba- 
tive calculation of the asymptotic amplitude. When the dipole mode decays, however, this 
term does contribute to the occupation of the residual state, with the asymptotic occupation 
probability of that state being 

Pdedb) = T / dt |ai||(M)| 2 (7) 



h 



— oo 
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b (he) 3 ' 

Using this perturbative probability, we can estimate the cross section for production of the 
decay product to be at least 

a dec » 2tt / Pdecip) bdb = 2irf3>y—^ (b max - b min ) , (8) 

Jb mln (he) 

where b min is the minimum value of the impact parameter at which the perturbative approx- 
imation is valid and b max = jv/uq is determined by limiting the impact parameter to values 
for which the adiabaticity parameter £ = ^ is less than onejfj, with hwo the minimum 
value of the excitation energy considered as contributing to the dipole mode. The resulting 
cross section thus grows with energy as 7 2 . 

The absurd result furnished by the above estimate is confirmed by the full coupled- 
channels calculations given in the table. The calculations were performed for the system 
208 Pb + 208 Pb using the three channels described above, with a dipole excitation energy of 
hw = 13.4 MeV, a reduced matrix element in accord with the giant dipole resonance sum 
rule and a minimum excitation energy of tvjjQ = 8 MeV. We have labeled the cross sections 
of the table as dipole-mode excitation cross sections. In the case of zero-width, they are in 
fact the cross sections for excitation of this mode. In the case of finite width, the values 
represent the flux that was excited to the dipole mode to later decay to the residual state. 

Similar calculations were reported in Ref. 2^|. Comparison of their zero- width results 
with ours and with those of Ref. |3| leads us to conclude that they performed the calcula- 
tions with the pure field form that we are advocating. The trend of their finite-width results 
also leads us to infer that they took the width into account in the decay of the dipole mode 
but not in its excitation. 
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Table I: Dipole mode excitation cross section of Pb incident on Pb using the pure field and 



the mixed (Lorentz) representations of the interaction. 



Representation 


T(MeV) 


cr(b) at 1 GeV/nucleon 


cr(b) at 10 GeV/nucleon 


field 





4.26 


13.98 


field 


4 


4.53 


13.73 


mixed 





4.46 


98.84 


mixed 


4 


15.09 


334.06 



We thus conclude that the spurious potential term in the interaction produces absurdly 
large cross sections when decay of the excited state is taken into account. Even when the 
width of the state is zero, this term, abetted by depletion of the ground-state, introduces 
contributions that increase with incident energy and also lead to absurdly large cross sections, 
such as those at 10 GeV/nucleon shown in the last column of the table. The unphysical 
results in both cases are attributable to the unphysical 1/6 term in the truncated coupled- 
channels calculations. We can also argue against such a term on simple physical grounds: 
we expect the polarization that produces the excitation of the target nucleus to depend on 
the gradient of the potential, that is the electric field, rather than the potential itself. In 
the following, we will show how the electromagnetic interaction Hamiltonian can be recast 
in a more physical form. 



Expansion of the interaction Hamiltonian 

We want to obtain the first few terms contributing to the interaction in the expansion of 
the electromagnetic fields about the center of the target nucleus, x = 0. Such an expansion 
is reasonable if the fields are slowly varying over the extent of the nucleus. 

We thus take 



j d 3 x p(x, t)ip(x, t) « J d 3 x p(x, t) ^po(t) + x ■ V(po(t) + ^xx ■ W(po(t) + 



(9) 



and 



J d 3 x J(x, t) ■ A(x, t) « J d 3 x J(x, t) ■ (A Q (t) + x ■ V 'A (t) + •••), (10) 
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where the subscript on the fields and their derivatives denotes the evaluation of these 
quantities at the point x — 0. 

Evaluation of the scalar potential terms is straightforward. Evaluation of the vector 
potential terms requires a bit more work. We use the continuity equation, 

V./+| = 0, (11) 



to obtain two supplementary identities [23]: 

J d 3 x J k = j d 3 x (V ■ (x k fj - x k V ■ fj = J d 3 xx k ^, (12) 

where the integral of the exact divergence is zero due to the finite extent of J, and, after a 
similar calculation, 

J d 3 x (JkXi + JiX k ) = J d 3 xx k x^. (13) 
Using the first of these, we can write 

/ d 3 x J(x,t)-A Q (t) = [ d 3 xJ k (x,t)A ok (t) = J2 j d 3 x^x k A 0k (t) = j d 3 x^x-A (t). 

J k J k J J 

(14) 

Using the second, we find, with a bit more work, 

j d 3 xJ(x,t)-(x-V)A (t) = \f d3x ^x- (x-V)A (t) + ^ J d 3 x (x x fj • (v x A (tjj . 

(15) 

Putting all the pieces together, we have 

V(t) = J d 3 x (p(x,t)ip(x,t) - ^J(x,t) ■ A(x,t)^j (16) 

= j d 3 xp(x,t) (ip {t)+x ■ Vip {t) + ~xx- VV</?o(*) + • • J 

-\ J d 3 x ^x ■ (l (t) + i(f • V)A (t) + ...) 

j d 3 x (fxj|^Vx ^o(*)) + • • • 

The unphysical long-range contribution to the excitation interaction of Ref. can be traced 
to the term containing the time derivative of the charge density, dp/dt. 
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The Gauge transformation 

We would now like to eliminate the terms containing the time derivative of the charge 
density, dp/dt, in Eq. (|TK|) . We can do this by making the gauge transformation A (of, t) given 
by the factor multiplying this term, 

A(x,t) = -x ■ ^A (t) + ^(x ■ V)A (t) + ...^ = -jT dux-A{ux,t). (17) 

The vector potential that results from the gauge transformation can be expanded as 

A'(x,t) = A(x,t)+VA(x,t) (18) 
= ~x x(vxl )+..., 

while the transformed scalar potential can be expanded as 

<f'{x,t) = <f(x,t)--— 19 
c at 

= I V ^(*) + --^ £ J +2 ff ' V ( V ^ o(t) + c^t ) + '"' 

We can then rewrite the interaction as, 

V(t) « y rf 3 x t) ^ (t) + ^ ■ fv^o(t) + J + • V ^V^o(t) + J + • 

A(fxj).(vxl (i))+.--, (20) 
which we can express in terms of the electromagnetic fields E and B at x = as 
V(t) « Jd 3 xp(x,t) \^p (t)-x-E (t)-^xx-VE (t) + --^ 



1 

2c 



1 

2^ 



<Tx [xx J) -B (t) + ■■■. (21) 



If we assume as well that the field-producing charge does not overlap with the nuclear one, 
we can write the electric quadrupole term as 

xx ■ VEq = XiXjdjE 0i = (xiXj — x 2 5ij/3)djE 0i , (22) 
since, in that case, V • En — 0. 
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We can then write the interaction in a form that parallels the classical (Cartesian) mul- 
tipole expansion |23], as 



V(t) = q<p(t) - d- E (t) - - QijdjE 0i (t) - m ■ B {t) + ■■■, (23) 
where q is the charge, 



2 

';3 



q = J d 3 xp(x ) t) 1 (24) 

d is the electric dipole operator, 

d = J d 3 xxp(x,t), (25) 
the Qij are the traceless electric quadrupole operators, 

Qij — J d 3 x (xiXj - x 2 5ij/3) p(x,t), (26) 
and rh is the magnetic dipole operator, 



rh = y J d 3 xx x J(x,t). (27) 

The multipole expansion does not depend on whether we have included magnetization cur- 
rents, exchange contributions or other corrections that determine the detailed structure of 
the target charge and current densities. The only property of these densities that we have 
used is the continuity equation, Eq. (fTT| . which should be valid in any case. 

The multipole expansion given here can be extended to all orders without great difficulty. 
It can be expressed compactly as 

V(t) = qip(t) — J d 3 xp(x,t)x- J duE(ux,t) — ~ J d 3 x |lx jj ■ J duuB(ux,t). (28) 

It satisfies the gauge condition x ■ A(x, t) =0, which can be rewritten as 

i 

dux-A(ux,t) = 0, (29) 

o 

(see Eq.[T2j) and is known as the multipolar or Poincare gauge j^. A drawback to our 
formulation is that it has been performed in a Cartesian rather than a spherical basis. 
However, the expansion in the spherical basis should be directly deducible from the Cartesian 
one, since the two are equivalent. 
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Discussion 



In the physically intuitive form given above, the electromagnetic interaction poses no 
problem, even when the excitation and decay of resonant states is considered. As an example, 
we take the Lienard-Wiechert potential due to a relativistic nucleus of charge Z passing on 
a straight-line trajectory with velocity v in the z direction at a distance 60 from the center 
of the charge/current distribution of the target. We then havej^J 

Z G< 1) 

ip(b, z,t) = 7 — , and A(b, z, t) = -z ip(b, z, t). (30) 

\J(b-b ) 2 + j 2 (z-vt) 2 c 

The transverse and longitudinal components of the electric field, which induce dipole exci- 
tation, are 

z & z & 

Eo±(t) = -jb - -^E and E o\\(t) = -ivt- (31) 

(b 2 + ( 7 ^) 2 ) (bl + (Tvt)*) 

and tend to zero as 6g 2 and &q 3 > respectively, as the impact parameter increases. At high 
energies, the excitation is dominated by the transverse modes, which produce a cross section 
that grows as hx(b max /Lj n ), as expected|f|. 

Bayman and Zardi|2J] have observed that the mixed representation of the interaction 
given in Ref. neglects relativistic corrections to the quadrupole and higher multipole 
terms, which become important at high energies. These are included and discussed in their 
work and in Ref. Q. They are also taken into account correctly (and automatically) when 
the pure field representation is used. 

As mentioned in the introduction, Bayman and Zardi have also recently compared calcu- 
lations of relativistic Coulomb excitation in the Lorentz and Coulomb gauges [Uj. They find 
that the cross sections in the Lorentz gauge increases dramatically with respect to those in 
the Coulomb gauge at energies above about 2 GeV/nucleon. Although the Coulomb gauge 
is not equivalent to the one we present here, we suspect that it too might include only phys- 
ical contributions to coupled-channels cross sections. We are of the opinion, however, that 
the optimal form of the interaction in such calculations is the one given here, in which the 
excitation interaction is expressed in terms of the physical fields and their derivatives. 
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Conclusion 

We have shown that, through an appropriate gauge transformation, the Coulomb excita- 
tion component of the semiclassical form of the electromagnetic interaction can be expressed 
in terms of the electromagnetic fields and their derivatives. Aside from making a satisfying 
parallel with the classical case, the pure field representation of the excitation interaction, 
known as the multipolar or Poincare gauge, provides physically reasonable cross sections, 
even when resonance excitation and decay is taken into account. 
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